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LOCAL EXISTENCE IN RETARDED TIME UNDER A WEAK DECAY 

ON COMPLETE NULL CONES 

JUNBIN LI AND XI-PING ZHU 


Abstract. In the previous paper |10| . for a characteristic problem with not necessarily 
small initial data given on a complete null cone decaying like that in the work [6] of the 
stability of Minkowski spacetime by Christodoulou and Klainerman, we proved the local 
existence in retarded time, which means the solution to the vacuum Einstein equations 
exists in a uniform future neighborhood, while the global existence in retarded time is 
the weak cosmic censorship conjecture. In this paper, we prove that the local existence 
in retarded time still holds when the data is assumed to decay slower, like that in Bieri’s 
work [5] on the extension to the stability of Minkowski spacetime. Such decay guarantees 
the existence of the limit of the Hawking mass on the initial null cone, when approaching 
to infinity, in an optimal way. 


1. Introduction 

In the previous paper m, we considered a characteristic problem of the vacuum Einstein 
equations, where part of the initial data is given on an asymptotically flat complete null 
cone. The decay rate given on the complete null cone in m inherits that in the remarkable 
works of the stability of Minkowski spacetime in [6] and [8], which means that the decay 
rate we considered, is the same as that of the complete null cones in the solutions in BM- 
We remark that, the decay rate in the above works does not necessarily satisfy the peeling 
properties, which are implied by a smooth conformal compactification, but a suitable notion 
of the future null infinity can still be defined. What we proved in m is, that if we start 
from a complete null cone with the above decay rate, then the solution to the vacuum 
Einstein equations will always contain a piece of the future null infinity. This gives the 
local existence in retarded time, and the global existence in retarded time is exactly the 
weak cosmic censorship conjecture, see the discussions in [3]. It was also proved by Cabet, 
Chrusciel and Wafo in [3], that if the complete null cone attaches the future null infinity 
smoothly (and therefore the peeling properties are satisfied), then the solution contains a 
piece of smooth future null infinity. 


The authors are partially supported by NSEC 11271377. The first author is also partially supported 
by the Fundamental Research Funds for the Central Universities. 
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In [2], Bieri extended the stability of Minkowski spacetime (due to Christodoulou and 
Klainerman 0), under weaker decay and regularity assumptions. Roughly speaking, 
Christodoulou and Klainerman [6] considered the asymptotically flat initial data {'S,g,k) 
with the following decay at infinity: 

gij = + 04(1’"^^^), hj = 03(r“®/^), 

where / = Om(r^) means d^f = for all k < m. Bieri considered the following: 

gij = 6ij + 03 (r"^/^), kij = 02(r"^/^). (1.1) 

Bieri’s result [2] requires one less derivatives and one less power of r. One of the common 
features of both proofs is making use of an optical function u, which increases towards 
future and is normalized at infinity, whose level sets Cu are outgoing null cones extending to 
infinity. Along these null cones, the geometry tends to being Minkowskian in suitable rate. 
For example, without rerard to the difference in the regularity considered in both works, 
the curvature component^, in the spacetime considered by Christodoulou and Klainerman, 
behave like: 

a, /3 = o(r“’^/^), p = 0 { r ~^), a = /? = 0(r“^rr^^^), a = 

( 1 . 2 ) 

and behave like the followings in the spacetime considered by Bieri: 

a, j3, p,a = (3 = 0{r~'^TZ^^‘^), a = 0{r~^TZ^^‘^)- (1.3) 

Here = 1 + As a consequence of less decay being assumed on the initial data, the 
decay of various geometric quantities induced on the null cones Cu becomes slower. 

The decay considered by Bieri is expected to be sharp. One reason is that, as pointed 
out in [2], borderline terms appear in the energy estimates, which suggests that any further 
relaxation on the decay may cause divergence of the integral and the bootstrap argument 
cannot be closed. We will make more comments at this point later in a clearer way. Another 
reason seems to be more reasonable, that is, Bartnik [T] proved that, the asymptotically 
flat initial data verifying (|l.ip always has a well-defined and unique ADM energy, and it 
is the optimal decay rate (without regard to the regularity) to guarantee the ADM energy 
to be uniquely defined (see HI)- 

We still give an additional reason. On a fixed complete null cone, the asymptotic be¬ 
havior (jl.3l) states as 

a,^,p,cr = o(r"®/2), /3 = 0(r"2), a = 0(r"^). 


^The exact asymptotic behaviors are written in an integral form. 
^See the next section for definitions. 
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The corresponding estimates for the connection coefficient^ state as: 

X = !]_ = o(r"3/2), 

Otrx — f^trx = 0(r“^), Otrx, f^trx = 0{r~^), x = 0{r~^). 


(1.4) 


Suppose that the null cone is foliated by spherical sections labeled by u, where r ^ u. 
Then we consider the Hawking mass of Su- 


m(u) = — 


1 + 




We can evaluate the derivative of m(u) with respect to u: 


dm r f 

du 1071 jg 


(Htrx - Htrx)^ - -Htrylxl^ + Htryj^p 




(1.5) 


where = — p + \{x^x) — di/kry is the mass aspect function relative to a transversal null 
direction orthogonal to Su, whose integral over Su determines the Hawking mass. By 
dm), the Hawking mass m(u) has a limit moo as ;u(and r) — >■ +oo. In addition, the decay 
conditions (jl.4p . especially the decay conditions on x and ry, guarantee the convergence of 
the Hawking mass in an optimal way. 


Therefore, it is natural to study the characteristic problem with the decay rates (USD, 
dm), but without assuming smallness. The main result of this paper is the following 
semi-global existence result (which is a rough form of Theorem [2|) : 

Theorem 1. We consider a characteristic problem, with the geometric quantities on the 
complete null cone behaving asymptotically like the following: 

X = o{r-^/^), C = o(r“3/2). 

X = 0(r"^), try-try = 0(r "2), try = 0(r"^); 

a,l3,p,a = o(r"^/2); 

/3 = 0(r“^), a = 0{r~^). 

Then the solution to the vacuum Einstein equations always exists in a uniform future 
neighborhood of the initial null cone. In particular, a global optical function u normalized 
at infinity is constructed. 

Remark 1. This theorem shows that, the weakest decay dOD, dm) can also propagate 
locally to the future in retarded time. Recall that, if the rescaled metrics on the spherical 
sections Su tend to a standard round metric along the null cone, then the limit of the 
Hawking mass is the Bondi energ'^. In the proof of the theorem, we do not need to assume 
that the spherical sections of the initial complete null cone Cq tend to a standard round 
metric, in other word, the limit of the Hawking mass is allowed to have nothing to do with 


^See the next section for the definitions. 

^See the discussions in m and the references therein. 
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the Bondi mass. In addition, after proving the semi-global existenee, we can also prove 
that limits of the Hawking mass on Cu, which is the level sets of the global optical function 
u, decrease as a function of u, using a similar argument of proving the Bondi energy loss 
formula in [8]. 


We end this section by making some comments of the proof. The basic idea follows from 
that of the previous work [10] . The major difference in the energy estimates is that the 
decay for a, (3, p and a is weaker than that in the previous work, we use a different group 
of weights to generate weighted energy estimates using the null Bianchi equations. For 
both pairs of the equations for IDa-Df} and I2.l3-D{p, a), we use the same weight which 
means that, we rewrite the null Bianchi equations as 

^(r^lapd^^) + D{2r‘^\l3\^dp^) = • • • , 

D{r‘^\l5\^dp^) + D{r‘^{\p\^ + \(j\^)dp^) = ■■■ . 

The right hand side consists of a precise divergence term on the spherical sections and 
lower order terms. By integrating over the whole spacetime manifold, we can obtain the 
desired estimates for the curvature components. Recall that we use the weight in |10] . 
Notice that this is a refinement of the multiplier vectorfield method used in [2], which relies 
on a carefully designed group of approximately Killing and conformal Killing vectorfields. 
For the angular derivatives of the curvature, we use the relation y ~ r~^ as before. 


As remarked before, some borderline terms appear in the energy estimates when we 
study the weak decay. Such terms appear in the energy estimates generated by the group 
I2.{p, (T)-Dj3 of the null Bianchi equations, that is, the terms with underline in the following: 

+ D{r‘^\§\^dp^) = ■ ■■ -r‘^{p{Ilx,a)+cr{^lx^^) . 

The norm of the underlined terms over the spacetime is estimated by (|| • || refers to 

II • Wlbs^,^)) 

2 

rp\\ ||(ry)*x||||a||du'du' 

i=0 

^ ru / j-e ru \ 1/2 

<^ / sup ||(ry dw'• ( / / Wrpfdudu'] 

Jo u \Jo Jo J 

The above three factors are expected to be bounded because the expected decay rates for 
X, p and a are x = o(r“^/^), p = and a = 0{r~^). We can observe that the 

norm of the underlined terms diverges if we do not use the full decay of x, p and a. This 
is exactly the reason why we call them borderline terms. 


1/2 

'^dtt' ] du' 


snp lyj ||a||^du'J 


\ 1/2 
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As the borderline terms appear in in the curvature estimates, similar phenomena appears 
in the estimate for the connection coefficients. Such terms appear in the estimates for 
and y^(ntrx), using the coupled system of the Raychaudhuri equation and null 
Codazzi equation: 

D{ntrx) = +‘2ui^tix, 

I II 

di^ i^x) = + f^x • 

'' Tfi ^ 

First of all, we need to integrate (11.611 from the last slice, therefore the convergence of the 
integral of the right hand side is essential. The coefficient — ^ before (fltrx)^ requires that 
both the terms I and II should behave like o{r~^) to guarantee the convergence. 

For the term /, this is true because it is exactly the way what x is expected to behave, 
which is one of the crucial conditions guaranteeing the convergence of the Hawking mass, 
see dn and m- Notice that the term I is exactly the term — ^fltrxIxP in dl-Sp . 
However, the estimate for (the third order angular derivatives) of x is derived using (II.7p . 
This derivation relies on /3’s full decay o(r~^/^), and r/’s full decay (the term HI) o(r“^/^), 
the latter of which happens to be another crucial condition responsible for the convergence 
of the Hawking mass. 

The behavior of the term II is determined by the behavior of a;, which measures the 
difference between the double null foliation and the geodesic foliation. In the current case, 
in order that II behaves like o(r“^), oj should behave at least like o(r~^). The top order 
angular derivative^ of ui should be estimated through the equation 

= • • • + d(/(f (fI/3), 

which only suggests that (the top order angular derivatives) of cu behaves like o(r“^/^), 
because /3 behaves only like o(r~^/^), in the current case. Apparently, this is not enough. 
However, the estimate for the next to the top order angular derivatives of u is derived 
through the equation 

D(jj = ***“(“ 

which suggests that uj behaves like o(r“^/^), since p behaves also like o(r“®/^), and it is 
assumed that the initial null cone Cq is foliated by geodesic foliation, that is, cu = 0. 
Similar losses in the estimates for top order quantities also appears in [5]. The above 
argument suggests that the top order angular derivatives of the connection coefficients 


we assume up to I order angular derivatives of the curvature components, then the top order angular 
derivatives of the connection coefficients refers to the order / + 1. 


( 1 . 6 ) 

(1.7) 
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may not be able to directly estimated purely under a double null foliatiorQ. Therefore, 
since no smallness is assumed, it is convenient to work under the assumption of the third 
order angular derivatives of the curvature component^. 

The remainder of this paper is organized as follows. Section [2] dehnes the double null 
foliation, the components of the connection and curvature under such foliation, and lists the 
null structure equations and the null Bianchi equations. Section [3] states the main theorem, 
which is Theorem [2l and outlines the main steps of the proof. Section U] is devoted to the a 
priori estimate. Theorem [3l which is the main step of the whole bootstrap argument. The 
Appendix sketches the proof of the construction of the canonical foliation on the last slice. 


2. Preliminary 


2.1. Basic Geometric Setup. We follow the geometric setup and notations in [5]. We 
use M to denote the underlying space-time (which will be the solution) and use g to denote 
the background 3+1 dimensional Lorentzian metric. We use V to denote the Levi-Civita 
connection of the metric g. 

Let u and u be two optical functions on M, that is 

gCVu, Vu) = g{Vu, Vu) = 0. 

The space-time M is foliated by the level sets of u and u respectively. Since the gradients 
of u and u are null, we call the the these two foliations together a double null foliation. 
We require the functions u and u increase towards the future. We use Cu to denote the 
outgoing null hypersurfaces which are the level sets of u and use C_^ to denote the incoming 
null hypersurfaces which are the level sets of u. We denote the intersection Su^u = Q.u'^Cu-, 
which is a space-like two-sphere. 

We define a positive function 12 by the formula 

12“^ = —2g{'Vu,Vu). 


^However, the Raychaudliuri equation ILI can be written in the following form: 

^trX = -i(trx')^ - Ixf 

where x' = S'Hd s is the affine parameter. The Raychaudhuri equation has a preferable form in 

geodesic foliation. This suggests that the top order angular derivatives of try can be estimated, but 
behaves differently between the original double null foliation and geodesic foliation, see [2] and also see [6]. 

^Under such a strong regularity, many other estimates are actually easier to derive, see m- It turns 
out that if we assume a strong regularity in the first place, many difficulties will not come out. 

®By appealing to the preferable form of the Raychaudhuri equation, combining the techniques in [9], 
which is also commented in m , we should be able to assume only the first order derivatives of the curvature. 
In this paper, we only focus on the problem of decay. 
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We then define the normalized null pair ( 63 , 64 ) by 

63 = — 2 nV;u, 64 = — 2 nVu, 

and define one another null pair 

L — ^^ 63 , L = 1264 . 

We remark that the flows generated by L and L preserve the double null foliation. On a 
given two sphere Su^u we choose a local orthonormal frame (ci, 62 ). We call (ci, 62 , 63 , 64 ) a 
null frame. As a convention, throughout the paper, we use capital Latin letters 
to denote an index from 1 to 2 , e.g. ba denotes either 61 or 62 . 

We define i;^ to be a tangential tensorfield if (/> is a priori a tensorfield defined on the 
space-time M and all the possible contractions of (j) with either 63 or 64 are zeros. We use 
D4> and Dsf to denote the projection to Su,u of usual Lie derivatives Cicf) and Cijf. The 
space-time metric g induces a Riemannian metric ^ on Su,u and f is the volume form of fj 
on Su,u- We use and to denote the exterior differential and covariant derivative (with 
respect to fj) on Su^u- 

We recall the definitions of null connection coefficients. Roughly speaking, the following 
quantities are Christoffel symbols of V according to the null frame ( 61 , 62 , 63 , 64 ): 

XAB = 5(Vyie4, es), r]A = -^5(V3eA, 64 ), to = ^f^ 5 (V 4 e 3 , 64 ), 

Xab = 9C^Ae3,eB), =-^9(V4eA, es), LI = ^f^fi'(V3e4, es)- 

They are all tangential tensorfields. We also define the following normalized quantities: 

x' = ^~^X, x' = C = ^(^ -^)- 

The trace of x and y will play an important role in Einstein field equations and they are 
denoted by 

i jAB 1 J.AB 

trX = $ XAB, try = § Xab' 

By definition, we can check directly the following useful identities : 

^logn = + T]), Dlog^} = Lo, ID log Q = u. 

We can also define the null components of the curvature tensor R: 

1 1 

OtAB B.( 6 a , 64 , Bb , 64 ), [3A ^ -^(^A; ^4 , ^3 5 ^4 ) ? P ^ ^ ’ 

—AB ^3? ^B T 3 ); ^A 2 ^ ^A, ^B 
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2.2. Equations. We first define several kinds of contraction of the tangential tensorfields, 
which are used in expressing the equations. For a symmetric tangential 2-tensorfield 9, 
we use 6 and tiO to denote the trace-free part and trace of 9 (with respect to ^). If 0 is 
trace-free, D9 and D9 refer to the trace-free part of D9 and D_9. Let be a tangential 
1-form. We define some products and operators for later use. For the products, we define 
{9i,92) = ^^^^^^{9i)ab{92)cd and (6,6) = ^'^^(6 )a(6)b- This also leads to the 
following norms \9\‘^ = {9,9) and = (66- We then define the contractions {9 ■ 6a = 

{9i ■ 92)AB = {(^i)a‘" {92)cb, 01 a 6 = {Qi)ab{(^2)cd and 6®6 = 6 ® 6 + 

6 ® 6 “ (6,6)^- The Hodge dual for ^ is defined by *^a = ^a^^c- For the operators, 
we define d(A^^ = ^^^a, cu/l^yi = ^^Cb and (di/(^0)yi = ^^9ab- We finally define a 

traceless operator (^g^AB = {fC)AB + {f^BA - dV^^^AS- 

The following is the null structure equations that are used in this paper, (where K is 
the Gauss curvature of 

L>(mrx) = -^(fitrx)^ - + 2a;0trx, 

RX = -a, 

DtiR = -io2(trx')^ - f^^lx'lg 
Dr] = n{x-V- P), 

Dfl = ^{x-V + P), 

Dui = n^{2{r],2)-\v\‘^ - p), 

Duj = n^{2{p,p) - \v\^ -p), 

1 

K = - -trxtrx + - (x, x) - P, 

d(^ (fix) = ^^i(^2trx) +nx-V- ^^trxv - ^{3, 

D{nx) = ^^{f + ^trXX - ^trxx), 

D(fltrx) = Xi^{2d\^p + 2\p\^ - (x,x) - ^trxtrx-f 2p), 

Dp = -n(x -9 + §) + 2^Iw. 

We also use the null frame to decompose the contracted second Bianchi identity V“RQ^..y 5 = 
0 into components. This leads the following null Bianchi equations!^ 

5a - ^fltrxa + 2wa -f fl{-y gg - (4r/ C,)®j3 + ^XP + ^X^^} = 0, 


^See Chapter 1 of [5] for the derivation of these equations. 
^See Proposition 1.2 of [5]. 
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Da — ^^tixa + 2uja + + (4^ — C)®^ + 3x/0 — 3*xc^} = 0, 

DI3 + ^r2trx/3 — Ox ■ fi — ojfi — 0{d^a + (?? + 2^) • a} = 0, 

D0 + -Otrx/3 — Ox ■ /3 — w/3 + 0{d^a + (?? — 2(^) • a} = 0, 

DJ3 + ^Otrx/3 - Ox • /3 + a;/3 - Oj^p + *^(T + Srjp + 3*r]a + 2x • ^} = 0, 

- Ox • ^ + + 0{^/9 - *^(T + Spp - 3*'qa - 2x • /?} = 0, 

Dp + ^Otrx/3 - 0{di/(^/3 + {2p + C, /3) - ^(x, «)} = 0, 

Dp+ ^Otrx/0 + 0{di/^^+ (27? - C,^) + ^(x,a)} = 0, 

Da + ^Otrxo- + Ojcu/l/? + {2p + (, */3) - ^x A a} = 0, 

Do- + ^Otrxo- + Ojcu/l^ + {2p - (, *P) + ^x A a} = 0. 

We denote the first order elliptic operators (or Hodge operators) appearing above in the 
null Bianchi equations, by Di,D 2 ; where 

Di : tangential one-form | i-A a pair of functions (d(A^^, cu/1^); 

D 2 : tangential symmetric trace-free (0,2) type tensorfield 6 1 —)• tangential one-form d(^0. 
It is easy to calculate the formal adjoint 

*Di : a pair of functions (/, g) i-A- tangential one-form —^/ -|- *^g] 

*D 2 : tangential one-form ^ 1 —)• tangential symmetric trace-free (0,2) type tensorfield 

We will denote any one of the above elliptic operators (or Hodge operators) and their 
formal adjoint by D,*D. 

Before proceeding further, we list the commutation formulas which are used for the 
estimates of derivative^]. 

Lemma 1. Given integer i and tangential tensorfield cf). we have 

i 

i=i 

^See Chapter 4 of [5] for the first group. The second group can be derived directly by the definition of 
curvature. 
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J = 1 
i 

i=i 

Here we use “■” to represent an arbitrary contraction with the coefficients by ^ or In 
addition, if (f) is a funetion, then when i = 1, all commutators above are zero; when i > 2, 
all i’s are replaced by i — I’s in above formulas. 


3. Main Theorem and Structure of the Proof 


3.1. The Statement of Main Theorem. We first introduce the following scale invariant 
norms (for 2 < p < oo and q = 1, 2): 


ll^llupfeu) = I / r 


n 


II'CI|]HI"(u,m) ^ II ) ^IIlp(h,m) ; 


Ij±1.3i2] 


I2il.’i2l 


i=0 

^ r“2 








In addition, we define 


IICIIl? XJfH" 

ljil.Ji2l 


/ sup r 


1 

q 
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The main theorem of this paper is the following. 

Theorem 2 (Main Theorem). Suppose that Cq and C_q are two interseeting null eones 
where Cq is outgoing eomplete and Cg is ineoming, and Sq = Cq f) Cg is a two-sphere. 
Suppose also that Cq and Cq are foliated by affine seetions whieh are labelled by two fune- 
tions s and s. Let A(s) and A(s) be the larger and smaller eigenvalue of r{s)~‘^^\s^Q with 
respect to and r{s) be the area radius defined by 47rr(s)^ = Area{Ss,o). If the 

initial data given on Cq |J Cq satisfy the following: 

Oq =sup|||rVx,r^/^£rtrx||H2(s,o) + lk^'^^CllH3(s,o) + trxllH3(s,o)} 

+ sup I ||r3/2£ r^/'^Chlu^s) } 

/ \ 

T* / T' \ 

+ sup A(s) + (inf A(s))“^ + sup --h ( inf --) < oo 

s ^ fil + S \ s I -\- s J 

7^g =||rV2a, 

+ Wr'^P, r'^P, r‘^CF, ra||L2]Hi3(o) 

+ sup ||r®/^/3, r‘^P\\m^s,o) 

S 

+ sup{||ra||H2(o,s) + lki2a||Hi(o,s) + \\rD^a\y{o,£)} < oo- 

S 

Then there exists an e > 0 depends on Oq, TZq and a global optieal funetion u such that the 
solution of vacuum Einstein equations exists in a global double null foliation 0 < n < +oo, 
0 < n < e. In addition, 0 —?> 1 as +oo. 

Remark 2. The last statement that n —>■ 1 as +cx) means that u tends to the affine 
parameter of the incoming null geodesic as +oo. This guarantees that the null cones 
Cu will not shrink when approaching to infinity. 

3.2. Structure of the Proof. We will prove the above Main Theorem in the remaining 
part of this paper. The structure of the proof is similar to that in m- We begin the proof 
by defining 

Me,A = {c > 0 : c satisfies the following two properties for small e > 0 and large A > 0}, 

where e is a small positive parameter and A is a positive large constant. They will be 
suitably chosen in the context depending only on Oo,IIq. 

(1) The solution of vacuum Einstein equations g exists in a double null foliation given 
by (n, n) for 0 < u < c, 0 < u < e, where u = s on Cg and u on defines a 
canonical foliation, which means the following equation is satisfied: 

logO = 0, f logf2 = + ^ Q((x,a) - (a,a)) - (p-p)^ , 


(3.1) 
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and tt = 0 on Sofl. 

(2) Written in the double null foliation given by {u,u), TZ < A, where 
K = .up { ||rV2„, 

+ sup {||r^/3, r^p, r^a, ra||L2H3(;u)} 

U 


We will prove, for e > 0 sufficiently small and A sufficiently large, = sup^e,A = +oo. 
The proof is divided into following steps: 

Step 1. In this step, the canonical foliation on Cq is constructed and therefore is not 

empty. Then we assume that = sup^ < +oo. 

Step 2. This step is devoted to the a priori estimate. It is not hard to see € As^a- We 

work on the space-time region Mu^^e which corresponds to 0 < u < u^, 0 < u < e. 

and then the function u restricted on C,, induces a canonical foliation. We will 

—* 

prove that, if e > 0 is sufficiently small, 

7^<C(Oo,7^o). 

In particular, we can choose A sufficiently large such that TZ <l A. 

Step 3. By the existence result in m, We extend the solution ^ to 0 < u < u^, + S, 
0 < u < £ + 5' for 6, 5' sufficiently small, and then construct a new optical function 
varying in [0,e], such that {u,us) is a new double null foliation, and us induces 
a canonical foliation on In addition, the norms TZ expressed in the new 

foliation are bounded by ^A. 

Therefore, we have u^ + 5 G As^a-, which contradicts to that = sup^^^A < +oo. Finally, 
a global optical function u can be constructed by a limiting argument. 

We will complete Step 2 in the next section. The proof of Step 1 and Step 3, which is 
the construction of the canonical foliation, will be sketched in the Appendix and the reader 
can refers to m and the references therein for the full details. 


4. The A Priori Estimate 


We define the following norms of the connection coefficients; 

C> = sup I r^llOtrx, V,!l\\u3{u,u) + r\\Z trx, trx, w||H3(n,«) | 

u,u ^ ^ 
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We will prove the following theorem, which is the content of Step 2. 

Theorem 3. Suppose that the assumptions in Theorem [H hold and the solution of the 
vacuum Einstein equations g exists in 0 < u < u^, 0 < u < e for some double null foliation 
given by two optical function u, u, and u induces a canonical foliation on C_^ . Then if e 
is sufficiently small depending on Oq and TZq, we have the estimates 

0,7^<C(Oo,7^o). 


The proof of the above theorem is divided into the following two propositions. 

Proposition 1. Suppose that the assumptions in Theorem\M hold and IZ < oo. Then for 
e sufficiently small, 

o<c{Oo,no,n). 

In particular, the estimates for r'^\\htTx,^\\u^{u,u) does not depend on 

IZ. 

Proposition 2. If e is sufficiently small, we have 

IZ<C{Oo,IZo). 


4.1. Proof of Proposition [T]. It suffices to prove the Proposition [T] under the following: 

Bootstrap Assumption: O < Aq 
where Aq is a sufficiently large number. 

Recall r is the area radius and let I{Su,u) be the isoperimetric constant of the spherical 
sections Su,u, we introduce an auxiliary bootstrap assumption: 

^0^\l + u)<r<AOo{l + u), ^ < 27r/(5,,„) < 40o. (4.1) 

4.1.1. Premilary Lemmas. Using the bound of the isoperimetric constant, we have the 
following Sobolev inequalities. 

Lemma 2 (Sobolev inequalities, see Section 5.2 of [5]). Given a tangential tensorfield f, 
we have for q € (2, +oo), 

II^ II */’lllH[l (m,m) ) 

II*/’IIl°°(;u,m) ^ II */’I|]H[2(m,m) ) 


We remark that the notation A < B means A < C{Oq,IZq)B. By the Sobolev inequal¬ 
ities, we can bound Irfltry;!, Irflxl < C{Oo). Then it is a direct consequence that, if e is 
sufficiently small, (14.ip holds with 20o in the role of Oq. Therefore (|4.1I) actually holds. 
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Using the bounds of isoperimetric constant and Gauss curvature we have the elliptic 
estimates. 

Lemma 3 (Elliptic estimates for Hodge systems, see Chapter 7 of [5]). Suppose that 

< C- 

Now assume that 6 is a tangential symmetric trace-free (0,2) type tensorfield with 

= f, 

where f is a tangential one-form. Then 

Assume that (p is a function with 
We have 

ll<(’llH3(n,u) ^ lk^/llHl(n,«) + II(n,u) • 

Here cp is the average of cp over Su,u ■ 

Remark 3. We remark we do not need to estimate ||</>||L2(n,M) first to apply the second part 
of the above lemma, which is not so in applying the first part. We avoid using the positivity 
of the Gauss curvature K. 


We also have the following Gronwall type estimates: 


Lemma 4 (Gronwall type estimates, see Chapter 4 of [5] or Chapter 4 of [8]). Suppose 
that there exists a constant c such that ^ c. Then for an s-covariant 

tengential tensorfield (p, 2 < q < +oo, and any real v, we have 


\r 4 ’\\'Li{u,u) ^ 


II^ (w* 


+ 


r 




(u/ 


and for e sufficiently small, we have 


ll<(’llL'?(u,?i) ^ Cq (||<(>||l'I(m,0) + II:D<(’IIl1L9(;u)) • 


In addition, we have the Gronwall type estimate for the derivatives along D direction; 
Lemma 5. For e sufficiently small, we have, for n = 2, 3, 

II<(’IIh"(u,«) ^ ll'^llH"(?i,0) + II:^'^IIl1H"(u)- 
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Proof. We apply the above Gronwall estimate on the following equation (for 1 < i < 3) 

i 

Df^(f = + [D, = f"D(i> + Y, 

i=i 

By Holder inequality and Sobolev inequalies, we have 

\\4>1 ■ 4>2\\Mi{u,u) ^ 

for i > 2. Therefore, 

ll'/’llH"(n,u) ^ ^ SUp ||nx||H"(«,M) SUp || . 

u u 

The conclusion follows by choosing e sufficiently small. I 


We remark that the above lemma holds with in the role of (j), r^^cf) in the role of 
^(f>, for any /r with the right hand side being finite, and then the constant depends on /r. 


We also have the following: 

Lemma 6. For e sufficiently small, we have, for n = 2, 

J \KD(f\\^2(^^y)dudu 

for any fi with the right hand side being finite. 


1/2 


Proof. By Lemma [5l we have 

Notice that the right hand side does not depend on u, we then take supremum with respect 
to u, and integrate the square of both sides over we obtain 


supllr'^ + 


/ ii2 

/ 

J U-^ u 

The second term in the right hand side is estimated by 


r-M- 




{u^u) 


du 


r (r \\r^-^/^mMr^i^,u)du) du<e r 

J \,d Q / J • 


r-M- 


=e 



r-M- 


0 J u. 




The conclusion then follows. 
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It is crucial here we can change the order of the double integration of the second term 
on the right hand side. But, in this paper, we only need to estimate 

\ 1/2 


, 1/2 






Notice that ||ri^(;/)||L 2 < ||ri^(?i>||L 2 therefore we have 

[3il.’i2l ' [“1.2.21 “ 






'io)+£Sup\\r>^D(t)\\l2 

,H2(u)* 


H ' Mil^r ,11... ,vy, - .. -..Ji^. . 

•±1’2L2] ^ ^ bil>U2] ^ ^ U hLl^-lL2\ 

We will also call the above two lemmas the Gronwall type estimates. 


We first establish the following lemma, which says that the geometric quantities share 
the same estimates up to a multiple by fl: 

Lemma 7. If e > 0 is sufficiently small, then 


||(ry)-MogO||Loc(„_„) + ||(ry)^logn||L 4 („_„) + ||(r’y)^logn||L 2 („^„) < r ^/'^C{Oo,TZo). 
In particular, we have C{Oo,TZo)~^ < — 1) < C{Oo,Ilo). 


Proof. Recall that on Su^^, H satisfies the equation 

logfl = 0 , f logfi = ^d]fp + ^ Q((X,X) - (X,X)) - (p-p)^ ■ 

The expression on the right hand side is invariant, in particular, does not depend on fl. 
Therefore, by elliptic estimate and Sobolev inequalities, we have 

||(ry)^ilogO||Loc(„^,o) + ll(r-y)'logII||L4(.,,o) + ||(ry f logQ||L2(,^,o) < r-^/%^^^,C{Oo,no). 

Because we do not change the foliation u on Cq, is extended as a constant along every 
null generator of Cq. The above estimates then hold along the whole Co, i.e., replacing 
by u for all 0 < u < . 


Then we can estimate || log 00 ( 5 ^ by the equation Dlogfl = w as follows: 

|logO| < |logn|s,ol + / k|dn'< oC(Oo,^0) + 

Jo 


and then multiply both sides by and choose e sufficiently small. The derivatives of 
log n can be estimated by commuting y/ three times to the equation D log 17 = w and then 
apply Lemma [5] to obtain the conclusion. I 


Remark 4. ITe make an important remark here. Because we construet the canonical 
foliation on the last slice C„ , 17 will not be constant 1 but only constant along the null 
generators of Cq. Therefore, the norms included in O, IZ written on u = 0 are different 
from the eorresponding norms in Oq, IZq. Fortunately, the above estimates for up to the 
third order derivatives of It ensure that the differenees are up to a constant depending only 
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on Oq, TZq, so we do not need to worry about this. In addition, the differences of some 
components do no exist, for example, fix (but not x itself), w, ij, p and a. 


4.1.2. Estimates for the connection coefficients. Now we turn to the estimates for the con¬ 
nection coefficients. For the structure equations along ID direction, we will apply Lemma 
[5] or Lemma [6] for f equals to up to the second or third order angular derivatives of the 
connection coefficients. 


As the first step, we consider the structure equations for Dy '. Diix' and Dj]. By 
commuting at most three times to those equations, we have 

||rx',rtrx',r^/%l|H3(n,n) ^C'(C’o) + Ik®,+ Ikx' • ^llLlH3(n) 

'^C{Oq) + r^^||L2]gi3(„) + eAq 

<cm 

if e is sufficiently small. 

By applying Lemma El we also have 

LS°H2(u) + e sup sup \\rr]\\j^2 H2(„) 

u^u U 

+ e sup ,H2(u) 

<C{Oo) + e sup I ||r^/^Ox||H 2 («,«) \\r^^'^r]\\u^u,u)} 

U,U ^ ^ 

+ CSUp||r3/2^||L2 H2(„) 

<C{Oo) 

if £ is sufficiently small. Here we use 

lk>llL 2 m 2 («) < sup ||r^+^/V||H 2 ( 2 i,«) 
which is derived by the Holder inequality. 


We then consider the structure equation 

Duj = n^{ 2 {p,p) - \p\‘^ -p) 

We also have 

Ik^a;||e3(„_„) ■ fl + r‘^\v\‘^\kiM3{u) + \\r‘^P\kiM^{u) 

<eAo + e^'^^l|T^p||L2H3(;u) k C{Oo). 

We also apply Lemma El to obtain 

<esup \\r^/‘^p ' 11 + + ^sup Hr^^VllLf^ ,h3(«) 
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<esup||r^r/-^ + r3|^P||H3(„_„) +esup||r^/VllL2 h 3(«) +e^o 

U,U U 

<C{Oo) 


We then consider the structure equations for Ujl a-iid the last of which 

is derived by the equation for ^(f^trx). These equations involve the first order derivatives 
of the connection coefficients on the right hand side. Therefore, we can only estimate at 
most the second order derivatives of ntrx and rj under the bootstrap assumption. By 
applying the Gronwall type estimates, if e is sufficiently small, we have 

lk^f^trx,rmrx||H2(n,„) <Oo + lk??||LiH3(«) + \\r'^\r]\‘^ 

+ lk^/^llLlH2(n) 

< 0 ( 00 ) + e sup ||rr?||H3(„,„) + eAg + ||rV||L2H2(«) 

U 

<cm 


•'^IIl1H2(;u) + I 11(m) 

<C{Oq,TIo) + eAq + §\\ijim^{u,u) 

<C'(Oo,7^o). 


We then use the structure equation for D_{^x) to estimate 
l|r'/'f^xllH2(n,«) <cm + e sup lk^'^NllH3(«,«) 

U 

+ e sup • ry, r^/^mryOx, r^/^mrxS^xIlHSfe^) 

U 

<C{Oo) + e sup ||r^/\||e3(„_„) + eAg 

U 

<cm 

if e is sufficiently small. We also estimate 

lk^^^^^XllL2 ,L“H2(n) ^<^(Oo) +eSUp||r^/\||L2 e3(„) 

+ 6 sup 11^3/2^ . T], ?’^'^^5^trxS^xllL2 ,H2(u) 

'll 

<C{Oq) + e sup \\rr]\y3(u,u) + ^ sup \\Prj ■ t], Pm^x^p r20trxf^xllH2( 

U^U u,u 

<C{Oo) 


if e is sufficiently small. 
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We then consider the structure equations along u direction. First of all, from the above 
estimates, we can actually choose e sufficiently small such that r^|Otrx| are bounded 

by a constant c depending on Oq, therefore the hrst part of Lemma U] is true. 


The equations along u direction should be integrated from the last slice . Therefore, 
we should first investigate the quantities on the last slice . Recall the equation which 
is satisfied by log O on the last slice: 

f log 0 = ^di^^ + i Q((x,x) - (£x)) - (P-P)) , 
which is coupled with log 0 = 0. Because ry + ry = 2^ log O, we have 

- (XW)) - (P-P) (4.3) 

on the last slice. This equation should be viewed as the equation for ry on the last slice. 


We can also deduce the equation satisfied by w on the last slice : 

2^^ = 2d)/(^(0/3) + d^(30x • ry + -Otr^ry) + OtrxdyA^ry + {F — F + Otrxp — Otrx • p) 

■(4.4) 


where p = p - ^(xw), 

F = ^Otrxp- + {(2ry - C,^) - + + ^trxIxP}, 

and w = —Otrx log 0 by log 0 = 0. 


To apply the elliptic estimates, we need the following lemma. 

Lemma 8. For e sufficiently small, 

< C’(7^o,Oo). 


Proof. The proof relies on the null Bianchi equations for Dfi and Dp. By the Gronwall 
type estimates, if e is sufficiently small, 

(pw),?’^^^x •^llLiH2(4i) 

<7eo + e^/2Ao7^ + e^/2||rV,rV||L2H3(„) < C(7eo), 

P\\w{u,u) ^'^0 + llr-^'^^Otrxp, r^/^ry • ffi ■ ffi r^/^x ' a||LiH2(«) + ||r-^^^^||LiH3(;u) 

<7^0 + + e^^^||r’^^||L2H3(u) ^ C{TZq), 

\\r'^§\\uHu,u) ^^0 + ||r-^(^^trX,^^X,w) • ' «IIlih2(;u) + lk«llLiH3(n) 
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<7^0 + e^'^^AoT^ + e^^^lka||L2H3(u) ^ C'(^o)- 


The estimate for K comes from the Gauss equation K + ^tr^trx 



= -P- 


We also estimate 

<7^o + e sup r liotrx, Ox, w, vWu^u^u) sup \\r^^‘^P, r^^'^cr, r^/^^||L 2 M 2 (n) 

U,U U IU,U*J 

+ esup||r3/2p,r3/V||L2 ]ei3(„) 

< 7 ^o + e(Ao + 7 ^) 7 ^ < C( 7 ^o) 


if e is sufficiently small. 


Now, we investigate Otrx on the last slice. The equation for ^(Otrx) on the last slice 
reduces to 

1 


^(Otrx) + ^OtrxOtrx = - 2p). 


We estimate 

„5^3 


r^y'*(Otrx)||L3(„^,,,) <Oo + Ijr^OtrxOtrx,r^|r/p + r^pO^||Lie3(„j 

<€>0 + eA^ + esup |||p||L2(ji^,n)||^ + !?llH2(n^,u)| ^ C{Oo,'JZo) 


We then consider the structure equation for I7(y^(Otrx)) which is obtained by com¬ 
muting with the equation for L)(Otrx). Notice that the commutator does not 

contain the third order derivatives of x a-ud trx, therefore we have 

||r®y^(Otrx)||L2(„,n) <C'(Oo,i^o)+ [ r2||r^y3(Otrx)||L2(;u',«)l|f^trx,a;||H2(„/,„)dM' 

J U_ 

+ [ »"^lk^y^(^^x)llL2(;u',«)l|f^X,^trx||H2(n',«)du' 

J u 


nu 

+ / 

J U 
ru. 


h 2 (u',m) + ll^trx||H 2 (n',M)||Ox, Otrx, w||]ei 2 („/^.u) 


/ r2||Otrx||L2(„',„)||r^y a;||L2(„,^„)du'. 


du' 


(4.5) 


We will firstly analyze the terms in the last two lines above. The third line is bounded by 
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The fourth line is crucial because f^tr^ itself will lose decay. Therefore the decay of 
should be good enough. By Holder inequality, we estimate 


Line 5 < 


Ju 


< 


1/2 


r-3/2c(Oo)du' ||r3y3^r5/2a;)||L2 l 2(.) < C(Oo). 


We now go to the equation for d^(Hx)- By elliptic estimate, we have 

||r3y^(Hx)||L2(«,„) < ||T^y^(5^trx)llL2(«,«) + Wr^x ■ V, rQtixri\\m^u,u) + \\rP\\m^u,u) (4-6) 
So, the second line is estimated as 

Line 2 < f l|r^y ^(mrx)||L2(„,„) ||f^£ 5^trx||H2(;u',«)dM' 

J U 


+ 


P- 
J u 


r ^/^||r^''^HxT^/%||H2(;u'«)+?’ ^\\r^i^xr^^‘^u}\^^u,u)+r ^||r'^'^^/3||H2(n,«) 


< 

r\j 


+ 


< 

rsj 


X ||r^/^Hx,r^/^Htrx||H2(n',n)du' 

i [ lk^y^(5^trx)||L2(„,„)||f^£5^trx||H2(«',«)dii' 

J U 

sup{||r^/2||Hx||H2(:u,«)lk^^NllH2(:u,«) + lkf^trx||H2(«,«) ,H2(n)} 

+ ,H 2 (n)l X sup llr^HtrxIl + ||r3/2Hx||L2 

■ / ||r^y^(Htrx)||L2(;u,«)l|f^£^^ti'XllH2(«',«)dM'+ C'(C’o,7eo). 

Ju 

Notice that we use ||r®/^;0||L2 h2('u') < lk^'^^;d||L 2 l°°h2 < C{TZq) by Lemma [8l Now, 

(|4.5I) becomes an integral inequality for ||r®y ^(Htrx)||L 2 (u,u)- To apply Gronwall inequality, 
we need the following: 

C—* _ ^ \ 1/2 

J ||H£Htrx,w||H2(„',«)dn'< (^y r-^^'^dvfj r^/^j^trx, r3/^a;||L2^ 

<C{Oo). 

Therefore, we have concluded that 

||r5y3(f^trx)||L2(.,.)<C(Oo,7^o), 
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and by Lemma [H (|4.2p and (14.61) . we have 

< C{Oo,n), 

||r 9 / 2 y 30 x||L 2 („) <C(Oo, 7 eo). 

Now we turn to ij. We apply the elliptic estimate to (14.31) on the last slice, 

\\r^^^v\\uHu„u) < \\r^^‘^p\\M^u,,n) + • xllH2(«„n) + Wr^^'^vh^u^u) < C'(C>o,7^o). 

We then apply the Gronwall type estimates to the eqnation for Notice that the 

commutator [D,'^^]r] = does not contain ^^rj, and then does not 

contain new terms that are not estimated before. Therefore, we have 

r^f^v\k^iu,u) 


P—* 

k\k^f^V\\TL‘2(u^,u) + f [II^XllH3(ti',ti)(lhllH2(;u',u) + lhllH3(u',n)) + II/3 ||h3(;u',-u)] dw' 

Jy^ 

P—* r 

\\r^x\\m3iu',u){\\r^^‘^v\\mHu',u) + lk^'^NllH3(«',«)) 

Ju ^ 

+ lk^^^^llH3(u',M) du' 

k\\r'^f^v\h{u„u)+r~^^‘^ sup|||rOx||H3(„,«)(lk^'^NllH2(«,«) + lk^'^\llH3(;u,«))} 

Therefore, multiplying both sides by we have 

lk"/"y\llL2(n,«) <C’(Oo,7^o,7^). 

Finally we turn to u. By elliptic estimate (Lemma [8] and [3]) , we have 

lkiLllH3(u,,«) ^ + \k^\\h'2{u^,u) k c{Oq,tZo). 

We then apply the Gronwall type estimates for z = 0,1, 2,3, we have 

/■“* r ^ — _ 

lllLllH3(M,n) ^||lil||H3(M^,u) + / II^X) ^trx||j|2('^/^„)||w||][|2('„/„) 

J U_ *- 

+ ||2(?],^) \q\ ||]H[3(^/^^) + l|p||Bl3(n',ti) dn . 

By the Gronwall inequality, the first term in the integral is absorbed and we then have 
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+ r h3(«)- 

Therefore, multiplying both sides by r, we have 

lk^llH3(n,u) ^ C{Oo,TZo,TZ). 

4.2. Proof of Proposition [2l The proof relies on the following energy estimates. 

(i) 

Lemma 9. There exist functions rj ' with z = 0,1, 2 and j = 0,1, 2,3 such that we have, 
for every u € [0,;u^], 

ll^^^"«llLf„_„,H3(n) + II^"/3 |Il2H3(.) „ ,H3(0) + V 






0 JO 


i=0 

ik^^^/3iiLf„^^]H3(«)++ ypy^whmm 



/o Jo ^^0 

Jo Jo ^^0 

+ WmWllMHu) ^ll^'^'/3|lLf,„,H3(0) + W^'^^WllMHO) 


^[0,u*] 




10 Jo 


with Tj ' satisfying the following estimates: 

3 

Y1 II'^o*^IIl 1(«,«) ^ll^^trX>X,w||||af + h,!Zllll“llll/5|l + ll^llH2(n,,.)I|a||ll^ll + llxll l|a|| l|p, 


cr 


2=1 


(4.7) 


^ lki*^llLi(«,«) <ll^^trX>X,w||||/3|P + ||r?,^||||/?||||p,o-|| + ||xlll|a||l|p, o-|| + ll^llH2(2i,t 

+ llxllll/3||||/3|| + ||mry,x||||p,af, 


a 


2=1 


(4.8) 
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i=l 

+ llxlllb,o-||||a|| + ||^^trx,£w""'^"^ 


+ \\K\\m^^u^u)\\p,o- 


(4.9) 


^ lk3*^llLi(«,«) <ll^^trx,X,^llll^f + ll??,^llll^lll|a|| + llxlll|p,f^lll|a|| + ll^llH^fe^)ll^lll|a 

i=l 

+ l|f^trx,X,w||||af, 

where the norm || • || above refers to || • ||h3(«,m)- 


(4.10) 


Proof. We can simply define the error terms Tq \ T 2 *\ Tg*^ to be 

r2+2v«d//^ =^(r2+2*|y*a|2d/i^) + D{2r^+'^^\f^ (411^• y*/3)d/r^, 


r2+2v«d^^ =^(r'+'iy*/3pdM^) + D{r‘^+^\\f^ p\^ + \f^a\^)<ip^) 

- 2r2+2M(A^(ll(y V . TP - T<y • y**/3))dii^, 
r2+2vfd;i^ =D(r2+2*(|y V|2 + ly V|2))d/i^) + Zl(r2+2*|y*^|2d^^) 

+ 2 r 2 + 2 M^(Q(y V • y */3 - y V • y“/ 3 ))d/i^, 

r^vfd/i^ =^(2r2‘|y*^|2d/r^) + Zl(r2‘|y*«|2d^^) + d(A^(41iy*a • y*^)dM^. 


(4.11) 

(4.12) 

(4.13) 

(4.14) 


The estimates (I4.7l) - (j4.10p then can be derived by commuting y* with the null Bianchi 
equations. Notice that the Gauss curvature K comes from the commutator of ^ and the 
Hodge operators. We omit the details here. Then the estimates ()4.7p - (j4.10p follows by 
integrating (14.IIP . (I4.12p . (|4.13p . (14.141) over Sw,u for all u £ [0,u^]. I 

(u^u)e[0,;u] X [0,£] 


( 2 ) (” 2 ) '--- 

One should notice that in the estimates of T 2 and Tg , Htrx appears instead of Htrx. 


By (14.7p . (14.Sp . we have 


3 1 


Jo Jo i= 0 fc =0 


< 

rsj 


sup ||Htrx,X,w||H3fe«) ^ 
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+ supr{\\rj,rj,x,x\\m3{u,u) + \\K\\M^u,u)) { / /3\\l2 

u,u \Jo [ 0 ,ii] ^ ^ 


\ 1/2 


+ sup ||mrx,£w||H3(«,„) / \\r^^‘^l3,r^/‘^p,r^^‘^a\\l2 ]H[3(n')‘^“' 


u,u 

<eon^ 


1/2 


By (14.9p . we have 

3 




0 ./o 


i=0 


sup ||S^trx||H3(„^„) ^ 

+ sup r{\\p, p, x\\m3 {u,u) + \\K \\m^ i^u,u)) ( [ h3(«')^“J 


X pe \ 1/2 

+ lk^^^XllL2,„jLg°H3(u)l|l’^'^^P)'l’^^^<7||L2^^]Hi3(„) ■ e^^^||ra||L2H3(„) 

+ supr3/2||Otrx,X,w||H3(„,«) / (1 + M')“^'^^lk^^llL2H3(«')dli' 

<e^l^On^ + 0^“(1 + n')-'/'l|r'^llL=H3(.odM' 

The factor (1 + in the last term above, which ensures the convergence, is exactly 

the reason why we need ntrx appears instead of ntrx in the estimate of . This is also 
the case in the estimate of Tg^, The term with underline which comes from the estimate 
of r 2 *^ is the borderline term. And this is exactly the place the norm comes in. 

By (I4.10p . we have 


^ r2+2* ||rW ||Li(«',«')dn'dM' 
i=0 

SUp||mrX,X,w||H3fen) [ „jH3(n')'^“' 


'0 x'O i=0 
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+ supr(||ry,^,x||H3(:u,«) + II^IIh2(«,«)) 

+ supr3/^||ritrx,X,w||H3(„,„) / (1 + lk«llL2H3(n')dM' 

u,u Jo u / 

(l+M')”^^^lk«llL2H3(nO^^^' 


\ 1/2 


\^l/2 


<e^l^OV? + O 


C^') C 2 ) 

By the first two estimates in Lemma[9]and the above estimates for Tq and t{ , we have 
I|i’^^^«>1’^/2/3|Il2h 3H + l|l’^/5>l’V,l’^l^llL2H3(n) ^ Tll^sOV?. 


By Proposition [H O < C{Oo-,'R-o-,'R)-, therefore if e is chosen sufficiently small, we have 

||r5/2a,r5/2^||2^jj3(„) + Wr"^^ '^(^o)- 

Remark 5. The differences between the initial norms ||r^/3, r^p, r^ci, r^/3, ra||L 2 jj 3 (Q) and 
the corresponding norms in TZq are determined by all third order derivatives (including ^, 
D and their mixed derivatives) of it, which are bounded by a constant depending only on 
Oq, TZo if £ is sufficiently small. 


By the last two estimates in Lemma [9] and the above estimates for T 2 *^ and 


3 ! 


I^^'^IIl2h 3(«) + ll™llL2H3(n) ^^0 +£^^‘^0TZ^ 

+ 0 f (1 + u') ^'^^(lk^^llL2H3(:u') + l|l’^llL2H3(«'))di^^- 


Notice that O appearing in the third term on the right hand side comes from the estimates 
for X) and w, which depend only on Oq, TZq. Therefore, by choosing e sufficiently 

small, and absorbing the last term, we have 


II^^^IIl2H3(;u) + ll’’^^llL2LjJ^jH3 ~ <^(^> 0 ,^ 0 )- 

(*) (0 

We then substitute this estimate back to the estimates for t) and Tg , and again by the 
last two estimates in Lemma [U we have 


IrVV, 




Therefore, we complete the proof of Proposition [2] and Theorem [3l 
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5. Appendix 

As mentioned above, we will sketch the proof of Step 1 and Step 3, which is the con¬ 
struction of the canonical foliation on an incoming null cone. We refer the readers to m 
and the references therein for the details. 

Recall that a canonical foliation on an incoming null cone means that, under this 
foliation, the following equation holds: 

logO = 0, f logf^ = + ^ Q((X,X) - iZx)) - (p-p)^ , 

which is exactly (13.11) . Now we work on an arbitrary incoming null cone with a back¬ 

ground foliation given by a function s. If we work on the initial slice Cq, the background 
foliation refers to that given by the affine function s, and if we work on the last slice 
Q-u -|-<S’ the background foliation refers to that given by u, which is obtained by extending 
the outgoing null cones to and gives a canonical foliation on , 

We use a function W = W{s, 6) defined on [0, e] x Su^ to represent a new foliation in the 
following way; The new foliation function is defined by the relation (s, 0),9) = 

s. Under the new foliation, we have also the new “lapse” function W represents a 

foliation iff 

os 

Here H is the lapse function relative to the background foliation given by s. In general, W 
does not necessarily represent a foliation but only a family of spherical sections parame¬ 
terized by s. 

Now given W > 0. We consider a map A, such that A{W) > 0 is again a family of 
spherical sections, which is dehned by the following: 

A{W){s, 9) = r 9),9)n-\W{s', 9),9)ds' (5.1) 

Jo 

where is the solution of the equation 

log 9), 9) = (^)G(iy(s, 9), 9), 

^'^^log = 0, 

where 

("'iG ^ i (^((‘"'’7. '"'k) - ‘“’'P'TTl)) - r'p - . 
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If VK is a fixed point of A, then the foliation given by W is canonical. The strategy is 
to hnd a suitable closed subspace of the functions IT, and prove that A restricted to this 
subspace is a contraction. 

We introduce the following norms on C_^ with respect to the background foliation given 
by s: 

0(:u) = sup ||rx,rtrx,r^/%,r^/%,rw||H 2 („_s) +sup||r^||]Hii(„_^) + sup 

s s s 

R{u) = sup + sup ||r^||ei(„^s) + sup ||r-^^a||L2(n,s)- 

s s s 

Using the argument in [TO], we can prove the following 

Proposition 3. Given two families of sections represented by ITi,IT 2 . T/0(;u),M(n) < C, 
then if e, Wi and sup^ ||r^/^(ry )^’^IT(s, •)||l 2 („^ 5 ), i = 1,2 are sufficiently small (depending 
on the initial data Oq, TZq and C), there exists a constant c depending on C such that 

||rV2(log ((^lU)O(ITi(s, •), •) - log •), •))llH2feO) 

< c\\r^/^{W,{s,-) -W2{sr))\\mHu,0)- 

Remark 6. Notice that the regularity of the norms in O and M are stronger than those 
used in HO). This is because we use three order derivatives of the curvature in this paper, 
therefore the second order derivatives of the curvature components are controlled. On the 
other hand, the function W — s, and the lapse — 1, decay only like (but not r~^) 

because p decays only like r~^l‘^ (but not r~^). 

The construction of the canonical foliation on Cg and is then a direct consequence. 

We discuss the more subtle case for and the case for Cq is similar. We introduce 

the closed subspace K, = C C{[0,s], H'^{Su^^s,o)) such that 

'W{0,9) = 0, 

< 0 <W{s,9) <s + 6', (5.3) 

SUPs ||r^/^(IT(s, •) — s)||]h[ 2(„ ^+5,0) < e/C, 

for some small number ex. to be hxed. Recall that s here is actually the origin optical 
function u, which induces a canonical foliation on Cg. The second condition ensures that 
the foliation given by IT € /C does not go beyond s € [0, e + d']. Notice that we can choose 
ex small enough such that the second condition holds. Then the existence of a function ug 
that induces a canonical foliation on follows from the following proposition; 

Proposition 4. Fore sufficiently small depending on Oq, TZq, and 5 sufficiently small (may 
depending on e and 6'), A{K,) C K, and A is a contraction in fC C. (^([O, e], 

The proof of the above proposition is similar as in nm, using m and (15.21) . provided 
that 0(;u^ + 6),'R{u^ + d) < C{Oo,TZ q). The estimates for 0,M are already done, except 
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a and its angular derivatives, Dm and its angular derivative, Da and its angular 

derivative, and I^a. The estimates for a are obtained in a similar way to p, as in Lemma[ 8 j 
The estimates for Da and its angular derivative can be obtained by commuting D and D 
with the null Bianchi equation for Da. After commuting derivatives, by taking into account 
the null structure equations and Bianchi equations, the right hand side contains at most 
third order angular derivatives of the connection coefficients and curvature components. 
Then we can apply the Gronwall type estimates. The estimate for D^a can then be done 
in a similar way, by commuting ^ with the equation for Da. The estimates for Dm and 
its angular derivatives, and ^m, rely on the equation (j4.4|) . which is written on . We 
commute ^ and D^ with this equation, and use elliptic estimate. We can then choose <5 
sufficiently small, to conclude that the estimates on still hold. 

Finally, using the argument as in m, if e is chosen sufficiently small, we can extend us 
back to the whole spacetime, in the way that the level sets of us are outgoing null cones, 
that are orthogonal to the spherical level sets of us on ^ 5 . 
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